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Abstract 

We present a new microscopic stochastic model for an ensemble of interact- 
ing investors that buy and sell stocks in discrete time steps via limit orders 
based on individual forecasts about the price of the stock. These orders de- 
termine the supply and demand fixing after each round (time step) the new 
price of the stock according to which the limited buy and sell orders are then 
executed and new forecasts are made. We show via numerical simulation of 
this model that the distribution of price differences obeys an exponentially 
truncated Levy-distribution with a self similarity exponent fx ~ 5. 
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I. INTRODUCTION 



In the last years a number of microscopic models for price fluctuations have been de- 
veloped by physicists Jl|-|6| and economists The purpose of these models is, in our 
view, not to make specific predictions about the future developments of the stock market 
(for instance with the intention to make a fortune) but to reproduce the universal statistical 
properties of liquid markets. 

Some of these properties are an exponentially truncated Levy-distribution for the price 
differences on short time scales (significantly less than one month) and a linear autocorrela- 
tion function of the prices which decays to zero within a few minutes . 



We present a new microscopic model with interacting investors in the spirit of |p|,|2|J14 
that speculate on price changes that are produced by themselves. The main features of the 
model are individual forecasts (or prognoses) for the stock price in the future, a very simple 
trading strategy to gain profit, limited orders for buying and selling stocks 0] and various 
versions of interaction among the investors during the stage of forecasting the future price 
of a stock. 

The paper is organized as follows: In section 2 we define our model, in section 3 we 
present the results of numerical simulations of this model including specific examples of 
the price fluctuations using different interactions among the investors, the autocorrelation 
function of the price differences and most importantly their distribution, which turn out 
to be (exponentially) truncated Levy distributions. Section 4 summarizes our findings and 
provides an outlook for further refinements of the model. 

II. THE MODEL 

The system consists of one single stock with actual price K(t) and N investors labeled 
by an index i = 1, . . . , N. In the most simplified version of the model the investors have 
identical features and are described at each time step by three variables: 
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Pi(t) The personal prognosis of investor % at time t about the price of the stock at time t + 1. 

Ci(t) The cash capital (real variable) of investor % at time t. 

Si(t) The number of shares (integer variable) of investor % at time t. 

The system at time t = is initialized with some appropriately generated initial values for 
Pi(t = 0), Ci(t = 0) and Si(t = 0), plus a particular price for the stock. 

The dynamics of the system evolves in discrete time steps t = 1,2,3,... and is defined 
as follows. Suppose time step t has been finished, i.e. the variables Kit), Pi(t), Ci(t) and 
Si(t) are known. Then the following consecutive procedures are executed. 
Make Prognosis 

Each investor sets up a new personal prognosis via 

P t (t + 1) = (xP(t) + (1 - x)K(t)) ■ e r % (1) 

where x G [0, 1] is a model dependent weighting factor (for the investor's old prognosis and 
the price of the stock) and are independent identically distributed random variables of 
mean zero and variance a that mimic a (supposedly) stochastic component in the individual 
prognosis (external influence, greed, fear, sentiments • ■ -, see also @). 
Make Orders 

Each investor gives his limit order on the basis of his old and his new prognosis: 
P(t + 1) - P(t) > 0: 

investor i puts a buy-order limited by Pi(t), which means that he wants to transform all 
cash d{t) into mt[d(t) / P^t)] shares if K(t + 1) < P(t). 
P(t + 1) - P(t) <0: 

investor % puts a sell-order limited by Pi(t), which means that he wants to transform all 
stocks into S^t) ■ K(t + 1) cash if K{t + 1) > P^t). 

Now let i±, i^i ■ ■ ■ , in a be the investors that have put a sell-order and their limits are Pi x {t) < 
Pi 2 (t) < ■ ■ ■ < Pi NA (t), and let ji, J2, • • • , Jn b be the investors that have put a buy-order and 
their limits are P h {t) > P h (t) >■■■> Pj Ng (t). 
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Calculate new price 

Define the supply and demand functions A(K) and B(K), respectively, via 

N A 

A(K) = J2S ia .0(K-P ia (t)) 

a=l 
N A 

B{K)=Y J ^S 3b -[l-6{K-P 3b {t))} (2) 

6=1 

with ASj b = mt[Cj h (t)/Pj b (t)] the number of shares demanded by investor j b , and 9(x) = 1 
for x > and 9(x) = for x < 0. Then the total turnover at price K would be 

Z{K) = min {A(K), B(K)} (3) 

and the new price is determined is such a way that Z(K) is maximized. Since Z(K) is 
a piece-wise constant function it is maximal in a whole interval, say K e [-Pi max , -Pj max ] for 
some ? max € {ii, . . . , in a } and j max G {ji, . . . , in b }- Then we define the new price to be the 
weighted mean 

+ 11 = ^OTb(70 ' (4) 

Note that the weight by the total supply and demand takes care of the price being slightly 
higher (lower) than the arithmetic mean (Pj max + -Pj max )/2 if the supply is smaller (larger) 
than the demand. 
Execute orders 

Finally the sell-orders of the investors ii,---,i max and the buy-orders of the investors 
ji, ■ ■ ■ , jmax are executed at the new price K(t + 1), i.e. the buyers ji, . . . , j max update 

S jb (t + l)=S jb (t)+mt[C jb (t)/P jb (t)} 

C jb (t + 1) = C jb (t) - K(t + 1) • (S jb (t + 1) - S jb (t)) (5) 
and the investors ii, . . . , i max sell all their shares at price K(t + 1): 

Si.it +1) = 

C ta (t+l) = C la (t) + S ta (t)-K(t+l) (6) 



If M p imJ < B(P jm J then investor j max cannot buy int[C jmaxW /P imax (t)] but only the 
remaining shares, whereas in the case A(P imax ) > P(Pj max ) investor i max cannot sell all his 
shares. The orders of the investors imax+i, ■ ■ ■ ,in a and j max +i, ■ ■ ■ ,3n b cannot be executed 
due to their limits. 

The execution of orders completes one round, measurements of observables can be made 
and then the next time step will be processed. 

A huge variety of interaction among the investors can be modeled, here we restrict 
ourselves to three different versions taking place at the level of the individual prognosis 
genesis: 

Ii: Each investor i knows the prognoses P^it), . . . ,P im (t) of m randomly selected (once 
at the beginning of the simulation) neighbors. When making an order, he modifies his 
strategy and puts in the case 

m 

m + 1) - [ 9i (t)m) + e 9i n (WM < (>)o (7) 

71=1 

a buy (sell) order limited still by his own prognosis Pi(t). We choose the weights 
gi(t) = 1/2 and g in (t) = I /2m for n = 1, . . . , m. 

I 2 : In addition to interaction Ii investor % changes the weights g after the calculation of 
the new price K(t + 1) according to the success of the prognoses: 

9i _(t + 1) = 9i _(t) - Ag 

g i+ (t + l)=g i+ (t) + Ag (8) 

where fro each investro i the index i_ denotes the investor from the set 

. . . , i m } with the worst (best) prognosis, i.e.: 

i- G {i, ii, . . . , i m } such that abs[P«_ (t) — K(t + 1)] is maximal 

i + G {i, ii, . . . , i m } such that abs[P i+ (t) — K(t + 1)] is minimal (9) 

The weight gi is forced to be positive, because an investor should believe in his own 
prognosis Pi(t). 
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I3: In addition to interaction I2 neighbors with weights gi_(t + 1) < are replaced by 
randomly selected new neighbors. 

III. RESULTS 

In this section we present the results of numerical simulations of the model described 
above. In what follows we consider a system with 1000 investors and build ensemble averages 
over 10000 independent samples (i.e. simulations) of the system. We checked that the results 
we are going to present below do not depend on the system size (the number of traders). 
When changing the system size, i.e. the number N of investors, the statistical properties of 
the price differences do not change qualitatively. Increasing N only decreases the average 
volatility (variance of the price changes). 

For concreteness we have chosen the following parameters: the initial price of the stock 
is Kq = 100 (arbitrary units, J7|), Each trader has initially Cj(t = 0) = 50000 units of cash 
and Si(t = 0) = 500 stocks (thus the total capital of each trader is initially 100000 units). 
The standard deviation of the Gaussian random variable z is o = 0.01 (with mean zero). 
We performed the simulations over 1000 time steps which is roughly 10 time longer than 
the transient time of the process for these parameters. In other words, we are looking at its 
stationary properties. 

First we should note that in the deterministic case a = no trade would take place [|TJ, 
hence the stochastic component in the individual forecasts is essential for any interesting 
time evolution of the stock market price. 

We focus on the time dependence of the price K(t), the price change A T (t) = K t+T — K T 
in an interval T, their time dependent autocorrelation 

r , , (A T (t + T)A T (t))(A T (t + T))(A T (t)) 

Ct(t> «a tW )s} - (A T ( t) r (10) 

and their probability distribution P(A T (t)). The statistical properties of the price changes 
produced by our model depend very sensitively on the parameter x in equation (TJ). In 
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particular for the case x — 1 it turns out that the total turnover decays like t~ 1//2 in the 
interaction-free case, which implies that after a long enough time no investor will buy or sell 
anything anymore. However, only an infinitesimal deviation from x — 1 leads to a saturation 
of the total turnover at some finite value and trading will never cease. 

In Fig. 1-4 we present the results of the interaction- less case with x = 1 (Fig. 1) and 
x = and contrast it with the results of the model with interactions ii, also for x — 1 (Fig. 
3) and x = (Fig. 4). 

For x = investor % does not look at his old prognosis but only at the actual stock price 
when making a new prognosis. In this case the distribution of the price can be fitted very 
well by a Gaussian distribution irrespective of the version of interaction or no interaction. 
The self similarity exponent 1/pL « 0.5 agrees with the scaling behavior of a Gaussian 
distribution. The autocorrelation function of the price differences decays alternating to zero 
within a few time steps. 

In the opposite case x — 1 investor i makes his new prognosis Pi(t + 1) based on his 
own old one and never looks at the current stock price. Now we can show that the distri- 
bution of the price differences decays exponentially in its asymptotic, but the self similarity 
exponent 1/// w 0.2 is too small to agree with a Levy stable distribution. The autocorrela- 
tion function of the price differences decays very quickly, so that there are significant linear 
anti-correlations only between consecutive differences. 





Io 


Ii 


h 


I 3 


x = 
x — 1 


0.442 
0.228 


0.466 
0.212 


0.472 
0.185 


0.472 
0.185 



The selfsimilarity exponent has been determined via the scaling relation P(A T = 0) ~ 
T -1 /^ and a linear fit to the data of P(A T = 0) versus T in a log-log plot. These least 
square fits yield the relative errors for our estimates of the self similarity exponent in 
the table above, which lay between 0.1% and 0.3%. 
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IV. SUMMARY AND OUTLOOK 



We presented a new microscopic model for liquid markets that produces an exponentially 
truncated Levy-distribution with a self similarity exponent ~ 0.2 for the price differences 
on short time scales. Studying the distribution on longer time scales we find that it converges 
to a Gaussian distribution. The autocorrelation function of the price changes decays to zero 
within a few time steps. The statistical properties of our prognosis oriented model depend 
very sensitively on the rules how the investors make their prognoses. 

There are many possible variations of our model that could be studied. It is plausible 
that a heterogeneous system of traders leads to stronger price fluctuations and thus a smaller 
value for the self similarity exponent /i (which appears to be ~ 0.7 for real stock price 
fluctuations fLOf). The starting wealth could be distributed with a potential law (comparable 
with the cluster size in the Cont-Bouchaud model). Or the investors could have different 
rules for making prognoses and following trading strategies. Another possible variation is 
to implement a threshold in the simple strategy in order to simulate risk aversion (the value 
of the threshold could depend on the actual volatility). 

Unfortunately, forecasts for real stock markets cannot be made with our model, because 
it is a stochastic model. We see possible applications for this model in the pricing and the 
risk measurement of complex financial derivatives. 

Acknowledgment We thank D. Stauffer for helpful discussions. H. R.'s work was 
supported by the Deutsche Forschungsgemeinschaft (DFG). 
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FIG. 1. Results of numerical simulations for the model without interactions Jo and x = (i.e. 
investors look only at their old prognosis Pi(t)). Shown are the price fluctuations for one sample 
(top), the autocorrelation function Ct(t) for T = 1 (middle) and the probability distribution 
P(At) of the price differences for T = 1. 
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FIG. 2. The same as Fig. 1, however with x = 1 (i.e. investors look only at the old price K{t)). 
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FIG. 3. The same as Fig. 1, however with interactions I\ (see text) and x = (i.e. investors 
look only at their old prognosis Pi(t)). 
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FIG. 4. The same as Fig. 1, however with interactions I\ (see text) and x = 1 (i.e. investors 
look only at the old price K(t)). Note the spikes in the time dependence of the price marking the 
significant enhancement of price fluctuations that lead to the truncated Levy-distribution of the 
price changes. 
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FIG. 5. The price fluctuations K(t) (top) and the price difference distribution P(A\) (bottom) 
of the model with interactions of the investors I2 (left) and I3 (right). The delta peak at Ai = 
comes from the events were no trade took place. 
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